Alternative Scenarios of Relativistic Heavy-Ion Collisions: I. Baryon Stopping 
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Simulations of relativistic heavy-ion collisions within the three-fluid model employing a purely 
hadronic equation of state (EoS) and two versions of the EoS involving deconfinement transition 
are presented. The latter are an EoS with the first-order phase transition and that with a smooth 
crossover transition. The model setup is described in detail. The analysis is performed in a wide 
range of incident energies 2.7 GeV < y^sjvjv < 39 GeV in terms of the center-of-mass energy. Results 
on proton and net-proton rapidity distributions are reported. Comparison with available data 
indicate certain preference of the crossover EoS. It is found that predictions within deconfinement- 
transition scenarios exhibit a "peak-dip-peak-dip" irregularity in the incident energy dependence of 
the form of the net-proton rapidity distributions in central collisions. This irregularity is a signal of 
deconfinement onset occurring in the hot and dense stage of the nuclear collision. 
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I. INTRODUCTION 

Over the last decade relativistic heavy ion physics has 
made tremendous progress in understanding of the QCD 
phase diagram in the domain of high temperatures and 
low baryon density However, a number of important 
questions still remain open. These are: At which incident 
energy an onset of deconfinement happen? What is the 
order of the deconfinement transition at high baryon den- 
sities? Is there a critical end point in the phase diagram? 
These questions form the main motivation for the cur- 
rently running beam-energy-scan program Q at the Rel- 
ativistic Hcavy-Ion Collider (RHIC) at Brookhavcn Na- 
tional Laboratory (BNL) and low-energy-scan program 
Q at Super Proton Synchrotron (SPS) of the European 
Organization for Nuclear Research (CERN), as well as 
newly constructed Facility for Antiproton and Ion Re- 
search (FAIR) in Darmstadt [1] and the Nuclotron-based 
Ion Collider Facility (NIC A) in Dubna @. 

This paper starts a series of papers, in which I hope 
to shed light on two first questions formulated above. 
In these papers I will report results of thorough simu- 
lations of relativistic heavy-ion collisions in the energy 
range from 2.7 GeV to 39 GeV (in terms of center- 
of-mass energy, y/s^lv). This domain covers the en- 
ergy range of the RHIC beam-energy-scan and SPS low- 
energy-scan programs, as well as energies of the future 
FAIR and NICA facilities and the Alternating Gradient 
Synchrotron (AGS) at BNL. Though experiments at the 
AGS have been already stopped, experimental data taken 
at the AGS are still unique since they were neither up- 
dated nor repeated in any newer measurements so far. 
The simulations were performed within a model of the 
three-fluid dynamics (3FD) @ employing three different 
equations of state (EoS): a purely hadronic EoS (hadr. 
EoS), which was used in the major part of the 3FD simu- 
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lations so far [8l-[ll| , and two versions of EoS involving 
the deconfinement transition [l^l . These two versions are 
an EoS with the first-order phase transition and that with 
a smooth crossover transition. Since neither of EoS's in- 
cludes a critical end point, these simulations do not touch 
the last question formulated above. 

This report is planned as a series of papers because 
it concerns a great number of bulk observables (rapidity 
and transverse spectra, flow observables and multiplici- 
ties) for various species and a large number of incident 
energies, their comparison with available data, and also 
illustrations of global evolution of collisions within differ- 
ent scenarios (i.e. EoS's). Analysis of whole the set of 
observables will be useful for revealing possible correla- 
tions in the energy evolution of these observables within 
different scenarios, which can be used as experimental 
indications of the deconfinement onset or its absence. 

It is reasonable to start this series with analysis of the 
baryon stopping because a degree of stopping of collid- 
ing nuclei is one of the basic characteristics of the colli- 
sion dynamics, which determines a part of the incident 
energy of colliding nuclei deposited into produced fire- 
ball and hence into production of secondary particles. 
The deposited energy in its turn determines the nature 
(hadronic or quark-gluonic) of the produced fireball and 
thereby its subsequent evolution. Therefore, a proper re- 
production of the baryon stopping is of prime importance 
for theoretical understanding of the dynamics of the nu- 
clear collisions. I will argue that certain irregularity in 
the incident-energy dependence of the baryon stopping 
may indicate an onset of deconfinement. 

The paper is organized as follows. In the next section 
im a brief survey of the 3FD model is presented: basic 
ideas and choice of parameters relevant for the present 
simulations. Properties of the EoS's used in the present 
simulations arc illustrated in sect. IIIII Predictions of pro- 
ton and net-proton rapidity distributions and their com- 
parison with available experimental data are presented 
in sect. IIVI Analysis of the form of these rapidity dis- 
tributions and its evolution with incident energy rise is 
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done in sect. |Vl In sect. I VII a summury of results is 
formulated. Results on baryon stopping in simulations 
with deconfinement transitions have been already briefly 
reported (without details) in Ref. [l3| . 



II. 3FD MODEL 

The 3FD model @ is a straightforward extension of 
the 2-fluid model with radiation of direct pions [T^ - [l6j 
and (2+l)-fluid model [13, The above models were 
extend in such a way that the created baryon-free fluid 
(which is called a "fireball" fluid, following the Frankfurt 
group) is treated on equal footing with the baryon-rich 
ones. A certain formation time r is allowed for the fire- 
ball fluid, during which the matter of the fluid propagates 
without interactions. The formation time is associated 
with a finite time of string formation. It is similarly in- 
corporated in kinetic transport models such as UrQMD 
[H and HSD 

Unlike the conventional hydrodynamics, where local 
instantaneous stopping of projectile and target matter is 
assumed, a specific feature of the 3FD is a finite stopping 
power resulting in a counter-streaming regime of leading 
baryon-rich matter. The basic idea of a 3-fluid approx- 
imation to heavy-ion collisions [13, HH is that at each 
space-time point x = (t, x) a generally nonequilibrium 
distribution function of baryon-rich matter, fhi {x,p), can 
be represented as a sum of two distinct contributions 



/br(a;,_p) = fpix,p) + ft{x,p), 



(1) 



initially associated with constituent nucleons of the pro- 
jectile (p) and target (t) nuclei. In addition, newly pro- 
duced particles, populating the mid-rapidity region, are 
associated with a fireball (f) fluid. Therefore, the 3-fluid 
approximation is a minimal way to simulate the flnite 
stopping power at high incident energies. 

The above assumptions are implemented into the for- 
mulation of the 3FD model as follows @. There is a 
set of continuity equations (reflecting the baryon num- 
ber conservation) 



dM^) = 0, 



(2) 



d,Tr{x) 



F;{x)+Frp{x), 
F^{x)+Fl(x), 
FLix)~Fr,ix) 

'(5^ (x — x' — Uf{x')t) 



(3) 
(4) 



(5) 



where the F^ are friction forces originating from inter- 
fluid interactions. F^ and F^ in Eqs. ©-(SI) describe 
energy-momentum loss of the baryon-rich fluids due to 
their mutual friction. A part of this loss \F^ — F^\ is 
transformed into thermal excitation of these fluids, while 
another part [F^ -\- F^) gives rise to particle production 
into the fireball fluid (see Eq. ([5])). F^^ and F^^ are asso- 
ciated with friction of the fireball fluid with the p- and 
t-fluids, respectively. Here r is the formation time, and 



u'^Ax') 



ul{x') + uX{x') 
|up(x') +ut(x')| 



(6) 



is a free-propagating 4-velocity of the produced fireball 
matter. Accordingly to Eq. (O, this matter gets formed 
only after the time span UpT upon the production, and 
in different space point x' — '{Jp[x') r, as compared to 
the production point x'. 

The nucleon-nucleon cross sections at high energies 
arc strongly forward-backward peaked. This fact, which 
originally served as justification for subdividing baryonic 
matter into target and projectile fluids, was used in (22j 
to estimate the friction forces, F^ and F^ , proceeding 
from only NN elastic scattering. Later these friction 
forces were calculated [11] based on (both elastic and in- 
elastic) experimental inclusive proton-proton cross sec- 
tions. In the present calculations the following form of 
the projectile-target friction is used 



K = ^^pIpI [« - <) Dp + {u'; + uX) Di 



(7) 



=p or t, p =t and t =p. Here, p| denotes a kind of 



"scalar density" of the p- and t-fluids (see below). 
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(8) 



where m^r is the nucleon mass 
Spt(spt"4m2,)]i/2/2TO 



rcl 



pt - my [u-p -h <) , 
"j^ is the mean relative ve- 



for a =p and t, where Jj^ = UaU^ is the baryon current 
defined in terms of baryon density Ua and hydrodynamic r 
4-velocity normalized = 1. Eq. © implies cr£;(spt)= / 

that there is no baryon-charge exchange between p- and ^' 
t-fluids, as well as that the baryon current of the fireball 
fluid is identically zero, =0. Equations of the energy- 
momentum exchange between fluids are formulated in 
terms of energy-momentum tensors T^'^ of the fluids 



locity ofthep-and t-fluids, and (Tp/E{s-pt) are determined 
in terms of nucleon-nucleon cross sections integrated with 
certain weights (see [13, HI, for details): 

o'p(spt)= / daNN~>NX ( 1 - cos e'en! ) , (9) 

"'ecm<7r/2 V Pin / 



da 



NN^NX 



ecm<7r/2 



1 



Si, 



(10) 



Here the integration is restricted to the forward hemi- 
sphere {0cm < T^/'^) of the ccnter-of-mass scattering an- 
gles 6'cm, Pin = (spt/4-m^)i/2 and Ei,, = Sp{^/2 arc the 
in-coming momentum and energy of the nucleon in the 
NN cm. frame, respectively, and Pout and i?out are the 
corresponding out-coming quantities. ap{spt) is nonzero 
at any physical Spt, as it is seen from Eq. At the 

same time, the (JEispt) quantity, which is responsible for 
the flreball production, is zero for Spt below the inelas- 
tic threshold. The overall i?^ factor in Eq. ^ controls 
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unification of p- and t-fluid into a single one, when their 
relative velocity gets small enough (for details see 0). 

The above friction ((T]) is a certain extension of that 
derived in [23| . The original derivation [23j was per- 
formed under assumption that baryon-rich fluids consist 
of only nucleons, and only proton-proton cross sections 
were used in ^ and (|10p . The extension is required 
because the original derivation [23| does not take into 
account: 

(i) various mesonic and baryonic species produced in the 
collision 

(ii) possible multiparticle interactions which are quite 
probable in the dense medium, 

(iii) possible medium modifications of cross sections and 
effective masses, and 

(iv) quark and gluon interactions, if dcconfinement oc- 
curs. 

In view of these uncertainties, it is reasonable to make 
provision for tuning the above friction. For this purpose, 
tuning factors £,{spt) in the scalar densities of the p- and 
t-fluids are introduced 



piispt) 



bar. 



Us 



pt) 



(11) 



where p^'"'', p™'^''', Pa and are scalar densities of all 
baryons, all mesons, quarks and gluons, respectively, de- 
fined in the conventional way. This quantities are sup- 
plied together with EoS. Factors like 2/3 and 1/3 in Eq. 
(fTTj) take into account the assumed scaling of cross sec- 
tions in accordance with the naive valence-quark count- 
ing. In view of above mentioned uncertainties of the es- 
timated friction, in Eq. different tuning factors are 
introduced for hadronic and quark-gluon phases: S^h and 
^q, respectively. 

The friction between baryon-rich fluids was fitted to 
reproduce the stopping power observed in proton rapid- 
ity distributions for each EoS. The results, together with 
those for formation time r and freeze-out energy density 
£frz, are summarized as follows: 



• Hadronic EoS with incompressibility K 
(hadr. EoS): 



C^(s) = l + 



In 



1/4 



190 MeV 



(12) 



r 2 fm/c, Efrz = 0.4 GcV/fm^ 



(13) 



The factor is not applicable here because of 
the pure hadronic nature of the EoS. The incom- 
pessibility K = 190 MeV is also chosen on the con- 
dition of the best reproduction of available data. 

EoS with the first-order dcconfinement transition 
[i3 (2-phase EoS): 



euis) = 1, eM 



60^. 



T = 0.17 fm/c, Efrz = 0.4 GeV/fm^ 



(14) 
(15) 



• EoS with crossover dcconfinement transition [l2| 
(crossover EoS): 

C^(,s) = l, e,^(.) = 200 ^, (16) 



r = 0.17 fm/c, Sirz = 0.4 GeV/fm^ 



(17) 



Within hadronic scenario (hadr. EoS) the friction has 
to be enhanced in order to reproduce the baryon stop- 
ping at high energies, Eiab > WA GeV. Though such a 
enhancement is admissible in view of above mentioned 
uncertainties, the value of the enhancement looks suspi- 
ciously high. Indeed, at ^s^vat = 17.3 GeV, i.e. at the 
top SPS energy, = 2.2. 

At scenarios with deconfinement transitions there is no 
need to modify the hadronic friction. This can be consid- 
ered as an indirect argument in favor of such scenarios. 
At the same time, the quark-gluon modification factor 
decreases with the energy rise, which is in agreement with 
our expectations that the quark-gluon friction should get 
weaker at high energies because of approaching to regime 
of the asymptotic freedom. 

Freeze-out was performed accordingly to the procedure 
described in @ and in more detail in p4 . [25| . The baryon 
stopping turns out to be only moderately sensitive to 
the freeze-out energy density Efrz ■ The freeze-out energy 
density was chosen mostly on the condition of the best 
reproduction of secondary particles yields. 

The formation time t also affects the baryon stopping, 
especially at high incident energies (top SPS and higher 
ones), when the fireball fiuid is well developed. In fact, t 
reduces the effect of the friction between the fireball and 
baryon-rich fluids. The larger r is, the later this friction 
starts to act and hence the weaker effect is produced by 
this friction. Therefore, the fitted value of t is essentially 
related to the strength of the fireball-baryon-rich friction. 

There are other full-scale (i.e. (3+l)-dimensional) ap- 
proaches to modeling nuclear collisions, which take into 
account the dcconfinement trasition. These have their 
advantages, as well as disadvantages as compared to the 
3FD model. The conventional hydrodynamical model of 
Rcfs. [26i . I27I docs such simulations in a very similar 
way but without taking into account incomplete stop- 
ping of colliding nuclei at the initial stage of the reaction. 
Therefore, such kind of simulations are justified only at 
moderately high energies. Another class of fluid models 
uses (hadronic) kinetic codes to "cook" the initial flreball 
which is subsequently considered within the hydro sim- 
ulation with possible deconfinement transitions [28l - l3l| . 
Such approaches disregard effects of dcconfinement tran- 
sitions at the stage of inter-penetration of colliding nu- 
clei, and hence cannot be used for analysis of the baryon 
stopping, which is the prime goal of this article. Kinetic 
models with a deconfinement transition, i.e. A Multi- 
Phase Transport model [3^ and a more consistent model 
of Parton-Hadron-String Dynamics [ssj . overcome this 
problem. Moreover, they avoid the problem of freeze- 
out inherent in hydrodynamic models. However, kinetic 
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models are able to treat only a crossover transition. The 
first-order transition remains beyond the scope of kinet- 
ics. Contrary to other hydrodynamic approaches, the 
3FD model can treat a deconfinement transition at the 
initial stage of the collision with due account of incom- 
plete stopping of colliding nuclei, though in essentially 
rougher approximation than that in kinetic models. At 
the same time the 3FD model is able to work with the 
first-order transition unlike kinetic models. 

In the present run of computations, higher incident 
energies were reached as compared with previous runs. 
This became possible because of implementation of an 
adaptive grid in the code. The size of the cell is made 
gradually larger with the expansion of the system pro- 
ceeded. Thus, when the system occupies the larger space, 
the code does not require a larger number of cells and, 
hence, a higher RAM memory. The adaptive grid docs 
not make the accuracy worse since spatial distributions 
become smoother at the expansion, that relaxes require- 
ments on the grid step. The adaptive grid made possible 
computations up to 62.4 GeV incident energy in terms of 
the cm. nucleon-nucleon energy, i.e. ^^snn- However, 
results for the top energy of 62.4 GeV are still not quite 
accurate, since an accurate computation requires unrea- 
sonably high memory and CPU time. This should be 
kept in mind when results for this energy arc displayed. 



III. EQUATIONS OF STATE 

Figure [T] illustrates differences between three consid- 
ered EoS's. The deconfinement transition makes a EoS 
softer at high densities. The 2-phase EoS is based on 
the Gibbs construction, taking into account simultaneous 
conservation baryon and strange charges. However, the 
displayed result looks very similar to the Maxwell con- 
struction, corresponding to conservation of only baryon 
charge, with the only difference that the plateau is 
slightly tilted, which is practically invisible. This invis- 
ible slope of the plateau results from plotting the pres- 
sure at the additional condition of strange density being 
equal to zero rather than at constant strange chemical po- 
tential. Application the Gibbs construction in hydrody- 
namical simulations silently assumes that the inter-phase 
equilibration in the mixed-phase region is faster than the 
hydrodynamical evolution. 

The 2-phase and crossover EoS's still differ even at 
very high densities. The latter means that the crossover 
transition constructed in Ref. [l2| is very smooth. The 
hadronic fraction survives up to very high densities. In 
particular, this is seen from Fig. [2] the fraction of the 
quark-gluon phase {Wqgp) reaches value of 0.5 only at 
very high energy densities. In this respect, this version 
of the crossover EoS certainly contradicts results of the 
lattice QCD calculations, where a fast crossover, at least 
at zero chemical potential, was found |34| . Therefore, 
a true EoS is somewhere in between the crossover and 
2-phase EoS's of Ref. 




FIG. 1: (Color online) Pressure scaled by the product of nor- 
mal nuclear density (no = 0.15 fm"^) and nucleon mass (mjv) 
versus baryon density scaled by the normal nuclear density 
for three considered equations of state. Results are presented 
for three different temperatures T = 10, 100 and 200 MeV 
(bottom-up for corresponding curves) . 



Figure [5] demonstrates that the onset of the decon- 
finement transition in the calculations happens at top- 
AGS-low-SPS energies. Similarly to that it has been 
done in (ssj , the figure displays dynamical trajectories of 
the matter in the central box placed around the origin 
r = (0,0,0) in the frame of equal velocities of collid- 
ing nuclei: |x| < 2 fm, \y\ < 2 fm and \z\ < "/cm 2 fm, 
where 7cm is Lorentz factor associated with the initial 
nuclear motion in the cm. frame. Initially, the collid- 
ing nuclei arc placed symmetrically with respect to the 
origin r = (0,0,0), z is the direction of the beam. At 
a given density ub, the zero-temperature compressional 
energy, e{nB,T = 0), provides a lower bound on the en- 
ergy density e, so the accessible region is correspondingly 
limited. In the case of the crossover EoS only the region 
of the mixed phase between Wqgp = 0.1 and Wqgp = 
0.5 is displayed, since in fact the mixed phase occupies 
the whole {s-hb) region. The e-hb representation is cho- 
sen because these densities are dynamical quantities and, 
therefore, are suitable to compare calculations with dif- 
ferent EoS's. 



Only expansion stages of the evolution are displayed, 
where the matter in the box is already thermally equi- 
librated, as a rule. The exceptions arc central collisions 
at y^SNN = 27 and 39 GeV, in which the matter in the 
box is not still thermalized in the beginning of the expan- 
sion stage. This non-equilibrium stage of the expansion 
is displayed by dashed lines in Fig. O The criterion of 
the thermalization is equality of longitudinal and trans- 
verse pressures in the box with the accuracy better than 
10%. Evolution proceeds from the top point of the trajec- 
tory downwards. Symbols mark the time intervals along 
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FIG. 2: (Color online) Dynamical trajectories of the matter 
in the central box of the colliding nuclei (4fmx4fmX7cm4fm), 
where 7cm is the Lorentz factor associated with the initial 
nuclear motion in the cm. frame, for central collisions of 
Au+Au at = 3.3, 4.9, 27 and 39 GeV (6 = 2 fm), 

and Pb+Pb at ^/sJTn = 6.4, 8.9 and 17.4 {b = 2.4 fm). The 
trajectories are plotted in terms of baryon density (ns) and 
the energy density minus ub multiplied by the nucleon mass 
(e — mNns)- Only expansion stages of the evolution are dis- 
played. Symbols on the trajectories indicate the time rate of 
the evolution: time span between marks is 1 fm/c. For the 
2-phase EoS (upper panel) the shadowed "mixed phase" re- 
gion is located between the borders, where the QGP phase 
start to raise {Wqgp ~ 0) and becomes completely formed 
{Wqgp = 1). For the crossover EoS (lower panel) the dis- 
played borders correspond to values of the QGP fraction 
Wqgp = 0.1, 0.5, 0.7 and 0.9. Inaccessible region is restricted 
by e{nB,T = 0) — mNns from above. 



the trajectory. Subtraction of the niNnB term is taken 
for the sake of suitable representation of the plot. The 
size of the box was chosen to be large enough that the 
amount of matter in it can be representative to conclude 
on the onset of dcconfinemcnt and to be small enough 
to consider the matter in it as a homogeneous medium. 
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FIG. 3: (Color online) The same as in Fig. [2] but for trajec- 
tories for all three different EoS's within the same frame. 



Nevertheless, the matter in the box still amounts to a 
minor part of the total matter of colliding nuclei. There- 
fore, only the minor part of the total matter undergoes 
the dcconfinemcnt transition at lOA GcV energy. 

As seen, the deconfinement transition starts at the top 
AGS energies in both cases. It gets practically completed 
at low SPS energies in the case of the 2-phase EoS. In the 
crossover scenario it lasts till very high incident energies. 

The trajectories for different EoS's are very similar at 
lower incident energies, as it seen from Fig. [3] At higher 
energies trajectories for dcconfinemcnt scenarios remain 
very similar, while the hadronic-EoS trajectories differ 
from those mentioned above and exhibit a peculiar be- 
havior. It happens because of a long (as compared with 
the interpenetration time of colliding nuclei) formation 
time of the fireball fluid, see (fTB]) . At the first stage the 
expansion proceeds when the fireball fluid has not been 
formed yet. Then the formation starts and the energy 
density (but not the baryon density) even slightly rises. 
When the formation is practically completed, the tra- 
jectory returns to its normal evolution - downward in 
energy and baryon densities. 



IV. PROTON AND NET-PROTON RAPIDITY 
DISTRIBUTIONS 

A direct measure of the baryon stopping is the net- 
baryon (i.e. baryons-minus-antibarions) rapidity distri- 
bution. However, since experimental information on 
neutrons is unavailable, we have to rely on net-proton 
(i.e. proton-minus-antiproton) data. Presently there 
exist experimental data on proton (or net-proton) ra- 
pidity spectra at AGS and SPS Ho-lil ener- 
gies. These data were analyzed within various models 
li i, m, [13, m m ISSI- The most extensive analy- 
sis has been done in (y, |47[ . Here I would like to repeat 
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FIG. 4: (Color online) Rapidity spectra of protons (for AGS energies, left block of panels) and net-protons (p — p) (for SPS 
energies, right block of panels) from collisions of Au+Au (AGS) and Pb+Pb (SPS) calculated within three considered scenarios. 
Experimental data are from collaborations E895 [H, E877 [la, E917 E866 and NA49 [i^-Iil. The percentage shows 
the fraction of the total reaction cross section, corresponding to experimental selection of events. For the E917 data [331 these 
are 0-5%, 5-12%, 12-23%, 23-39% and 39-81%. Feedback of weak decays into p and p yields is disregarded. 



this analysis. The motivation is to perform simulations 
with different EoS's within the same dynamical model, 
i.e. the 3FD model, in order to reveal differences pro- 
duced by different scenarios. 

Figure |4] presents calculated rapidity distributions of 
protons (for AGS energies) and net-protons (for SPS en- 
ergies) and their comparison with available data. Notice 
that difference between protons and net-protons is negli- 
gible at the AGS energies. At the top AGS energy of lOA 
GeV their difference is 0.03% at the midrapidity, see com- 
pilation of experimental data in Ref. (soj . Gontribution 
of weak decays of strange hyperons into proton yield was 
disregarded in accordance with mesurement conditions of 
the NA49 collaboration. At the AGS energies the con- 



tribution of weak decays is negligible. Correspondence 
between the fraction of the total reaction cross section 
related to a data set and a mean value of the impact 
parameter was read off from the paper [FH in case of 
NA49 data. For Au+Au collisions it was approximately 
estimated proceeding from geometrical considerations. 

As seen from Fig. |4l at lower AGS energies all EoS's 
predict the same results, since at these energies only 
hadronic parts of all EoS's are relevant. Results of the 
2-phase EoS start to differ from those of the hadr. and 
crossover EoS's beginning from 6 A GeV and first in cen- 
tral collisions. At higher energies this difference extends 
to more peripheral collisions. Unlike other scenarios, the 
2-ph.-EoS distributions exhibit a dip at midrapidity even 
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FIG. 5: (Color online) Rapidity spectra of protons (upper panel) and net-protons (lower panel) from central collisions of Au+Au 
{b = 2 fm) at low RHIC energies. Feedback of weak decays into p and p yields is taken into account. 



in central collisions. This dip contradicts the available 
experimental data and is very robust: variation of the 
model parameters and ([TS]) in a wide range does not 
remove this dip. Therefore, it is a direct consequence of 
the onset of the first-order phase transition, which starts 
precisely at these energies in the 2-phase scenario, see 
Fig. [21 Calculations within one-fluid (i.e. conventional 
hydrodynamics) p6l [27| confirm this conclusion. This 
dip survives even in one-fluid calculations involving the 
Ist-order phase transition in spite of immediate baryon 
stopping inherent in the one-fluid model. Notice that 
one-fluid calculations without deconfinement transition 
manifest a "normal" (for the one-fluid model) result, i.e. 
no dip. Accordingly to analysis of Ref. [1^ this dip is a 
consequence of larger pressure gradients in the longitudi- 
nal direction developed in the deconfinement-transition 
scenario. In 3FD calculations, this dip transforms into 
midrapidity peak at higher energies (30A GeV and 4QA 
GeV). With further energy rise {Eiab > 40^ GeV) the 
midrapidity peak again turns into a dip, see also Fig. [S] 
The latter dip is already a normal behavior which results 
from incomplete stopping of baryons and takes place at 
arbitrary high energies. 

As has been already mentioned, the behavior "peak- 



dip-peak-dip" in central collisions within the 2-phase- EoS 
scenario is very robust with respect to variation of the 
model parameters and (fT5)) in a wide range. It cer- 
tainly disagrees with data at 8 A GeV, 10 A GeV and 40 A 
GeV energies. It also disagrees with data at 20A GeV 
and 30^ GcV. However, the latter data have preliminary 
status, and hence it is too early to draw any conclusions 
from comparison with them. This behavior is in con- 
trast with that for the hadronic-EoS scenario, where the 
form of distribution in central collisions gradually evolves 
from peak at the midrspidity to a dip. The case of the 
crossover EoS is intermediate. One could conclude in fa- 
vor of a weak wiggle, since the distributions at lOA and 
15 A GeV exhibit a shallow dip while at 20 A GeV looks 
like a plateau. 

Beginning from 158j4 GcV to higher incident energies 
(see Fig. [5]) predictions of different scenarios for the net- 
proton distributions remain quite similar, at the expense 
of the substantial enhancemeirt of the hadronic friction in 
the case of the hadronic EoS, see Eq. p^ . At the same 
time difference of proton spectra increase with the energy 
rise. In calculations for energies above 158A GeV contri- 
bution of weak decays into proton and net-proton yields 
were taken in to account in accordance of experimental 
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FIG. 6: (Color online) Midrapidity value of the net-proton 
(upper panel) and proton (lower panel) rapidity spectrum 
from central Au+Au (at AGS and RHIC energies, b = 2 fm) 
and Pb+Pb (at SPS energies, b = 2.4 fm) collisions as a func- 
tion of the incident center-of-mass energy for three considered 
EoS's. Experimental data are from collaborations E895 3hI. 
E877 [33l, E917 [H, E866 [H, NA49 [UH] and STAR [HI. 



procedure of STAR and PHENIX collaborations. 

Comparison with available data indicate certain prefer- 
ence of the crossover EoS, though the crossover scenario 
does not perfectly reproduce the data either. Predic- 
tions of different scenarios for net-protons diverge to the 
largest extent in the energy region 8^ GeV < Eiab < 
40A GeV. Unfortunately data at 20A and 30^ GeV still 
have preliminary status and disagree with any considered 



scenario. Updated experimental results at energies 20^4 
and 30A GeV are badly needed to pin down the prefer- 
able EoS and to check a trend of the "peak-dip-peak-dip" 
irregularity in the net-proton rapidity distributions. 

Preference of the deconfinement-transition scenarios is 
seen at incident energies above the top SPS one, see Fig. 
ini where midrapidity values of the net-proton and pro- 
ton rapidity spectra from central collisions of Au-fAu (at 
AGS and RHIC energies) and Pb-^Pb (at SPS energies) 
are plotted as functions of the incident center-of-mass 
energy. In Fig. |6] the midrapidity values are displayed 
in a wider energy range. I would like to remind that 
results for top calculated energy of ^/sNN = 62.4 GeV 
are very approximate, since a more accurate computation 
requires unreasonably high memory and CPU time. As 
seen, a visible difference between net-protons and pro- 
ton data, as well as between predictions of hadronic EoS 
and EoS's with deconfinement transitions starts only at 
RHIC energies. At these energies the hadronic scenario 
certainly overestimates the proton midrapidity density. 



V. ANALYSIS OF "PEAK-DIP-PEAK-DIP" 
IRREGULARITY 

Preliminary results of the above-discussed "peak-dip- 
peak-dig" irregularity have been already reported in 
Refs. m, [S^l. There the friction forces for the 2-phase 
and crossover scenarios were poorly tuned and hence the 
corresponding simulations poorly reproduced available 
experimental data. Therefore, conclusions were based 
on a certain trend of the results of simulations. Here I 
present calculations with thoroughly tuned friction forces 
in the quark-gluon phase, which made it possible to rea- 
sonably (and often better than in the hadronic scenario) 
reproduce agreat number of observables in a wider (than 
before [11,113]) incident energy range 2.7 GeV < ^/snn < 
39 GeV. 

In order to quantify the "peak-dip-peak-dip" irregu- 
larity, it is useful to make use of the method proposed 
in Refs. [ll, [H, [13]. For this purpose the data on net- 
proton rapidity distributions are fitted by a simple for- 
mula 



dN / r / / ^ / 

— = a(exp{-(l/ws)cosh(?; 



Vc 



+ cxp{-{l/ws)coah{y-'ycr,i + ys)}) (18) 



where a, j/s and Wg are parameters of the fit. The form 
(|18p is a sum of two thermal sources shifted by zti/s from 
the midrapidity. The width Wg of the sources can be 
interpreted as Wg ~ (temperature) /(transverse mass), if 
we assume that collective velocities in the sources have 
no spread with respect to the source rapidities ij/g. 

The above fit has been done by the least-squares 
method. Data were fitted in the rapidity range \y — 
ycm\/Vcm < 0.7. The choice of this range is dictated by 
the data. As a rule, the data are available in this rapid- 
ity range, sometimes the data range is even more narrow 
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FIG. 7: (Color online) Parameters of fit (|18|) . j/s (upper raw of panels) and Ws (lower raw of panels), deduced from analysis 
of experimental data of net-proton rapidity distributions (points with error bars), as well as from that of results of 3FD 
simulations with different EoS's (shaded areas between respective lines). Fits of preliminary experimental data are displayed by 
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Grey bands in all the panels indicate the areas: between 0.25 ycm and 0.35 y cm for Ws (lower raw of panels), and OAycm and 
0.5 ycm for j/s (upper raw of panels). 



(80 A GeV and new data at 158 A GeV The above 

constraint is imposed in order to treat different data in 
approximately the same rapidity range. Another reason 
for this cut is that the rapidity range should not be too 
wide in order to exclude contribution of could specta- 
tors. I keep the old data at 158A GeV H^l in the anal- 
ysis because these are known in a wider rapidity range 
as compared with the new ones [4^ . A narrow rapidity 
range results in large error bars of the fit. To evaluate 
errors of the fit parameters, I estimated the errors pro- 
duced by the least-squares method, as well as performed 
fits in different the rapidity ranges: \y — ycm\/ycm < 0.5 
and \y — ycm\/ycm < 0.7, where it is appropriate. The 
error bars present largest uncertainties among mentioned 
above. 

Similar fit was applied to the calculated distributions. 
Since experimental data at AGS and RHIC energies were 
taken from Au-|-Au collisions while at SPS the Pb-|-Pb 
collisions were studied, the calculations were performed 
respectively for Au-|-Au {b ~ 2 fm) and Pb-|-Pb (6 = 
2.4 fm) central collisions. In fact, at the same inci- 



dent energy the computed results for Pb+Pb collisions 
at 6 = 2.4 fm are very close to those for Au-f-Au at 6 = 2 
fm. Therefore, the related irregularity of the energy de- 
pendence of the fit parameters is negligible. Similarly to 
the experimental data, the fit of the computed results 
was also performed in two ranges, \y — ycm\lycm < 0.7 
and \y — ycm\/ycm < 0.5, in order to estimate uncertainty 
associated with varition of this range. This uncertainty 
turned out to be a dominant one in the case of computed 
data. Therefore, in Figs. [7]and[5]results of the fit of com- 
puted spectra are presented by shaded areas with borders 
corresponding to the fit ranges \y — ycm\/ycm < 0.7 and 

\y - ycm\/ycm < 0.5. 

Parameters ys and Wg deduced from the fit of experi- 
mental data exhibits no significant irregularities in their 
energy dependence: they monotonously rise with the en- 
ergy within the error bars of the fit, see Fig. [T] Grey 
bands in all the panels of Fig. [7] are drawn to guide 
an eye. They indicate the areas: between 0.25 j/cm and 
0.35 ycm for Wg (lower raw of panels) , and between 0.4 ycm 
and 0.5 ?/cm for ys (upper raw of panels), where the ma- 
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FIG. 8: (Color online) Midrapidity reduced curvature [see. Eq. (|19p ] of the net-proton rapidity spectrum as a function of the 
incident center-of-mass energy of colliding nuclei as deduced from experimental data and predicted by 3FD calculations with 
different EoS's. Upper borders of the shaded areas correspond to fits confined in the region of \y — ycm\/ycm < 0.7, lower 
borders - \y - ycm\/ycm < 0.5. 



jor part of experimental points are located. In particular, 
similar (but based on different, double-gaussian fit) anal- 
ysis of Ref. [IHl also shows absence of any spectacular 
irregularities in exitation functions of these parameters. 
This is more so in view of the fact that the analysis of 
Ref. [ii] was performed only at Eiab > 20A GeV. The 
parameters deduced from the fit of distributions com- 
puted within hadronic and crossover scenarios also man- 
ifest quite monotonous behavior with incident energy. At 
the same time, the results of 2-phase scenario exhibit cer- 
tain, however not strong, irregularity. 

The representation in terms of ys and Ws is not quite 
spectacular. These parameters are interrelated to some 
extent. They produce a similar effect on the rapidity 
distribution, especially if it is fitted in a narrow rapidity 
range. Therefore, it is desirable to find a single quantity 
which characterizes the shape of the rapidity distribution 
in the midrapidity range. Such a parameter is a reduced 
curvature of the spectrum in the midrapidity defined as 
follows 




= iUcm/wsf (sinh^ ys ~ Ws coshy^) . (19) 

The factor 1/ {ycmdN/ dy)^^^ is introduced in order to 
get rid of overall normalization of the spectrum, i.e. of 
the a parameter in terms of fit (fT8|) . The second part 
of Eq. ()19|) presents this curvature in terms of param- 
eters of fit (fTS]) . Thus, the reduced curvature, Cy, and 
the midrapidity value of the distribution are two inde- 
pendent quantities quantifying the the spectrum in the 
midrapidity range. Excitation functions of Cy deduced 
both from experimental data and results computed with 
different EoS's are displayed in Fig. HI Fig. [5] demon- 



strates how the distribution shape evolves from a convex 
form at low incident energies to a concave form at high 
energies. 

The irregularity in data is distinctly seen here as a 
zigzag irregularity in the energy dependence of Cy. Of 
course, this is only a hint to irregularity since this zigzag 
is formed only due to preliminary data of the NA49 col- 
laboration. A remarkable observation is that the Cy 
energy dependence in the first-order-transition scenario 
manifests qualitatively the same zigzag irregularity (mid- 
dle panel of Fig. as that in the data fit, while 
the hadronic scenario produces purely monotonous be- 
haviour. The crossover EoS represents a very smooth 
transition, as mentioned above. Therefore, it is not sur- 
prising that it produces only a weak wiggle in Cy. 

As it was explained in detail in Ref. [l^l, the "peak- 
dip-peak-dip" irregularity is very natural in a system un- 
dergoing a phase or crossover transition. First, it is as- 
sociated with the softest point of a EoS [HBl. Therefore, 
the irregularity is weaker in the crossover scenario than in 
the first-ordcr-transition one. Indeed, the softest points 
in the crossover EoS is less pronounced than in the first- 
order-transition one [s?! . There is no softest point in the 
hadronic EoS and hence there is no irregularity. 

The second reason of this irregularity is a change in 
the nonequilibrium regime. The 3FD model takes into 
account the leading nonequilibrium of the nuclear colli- 
sion associated with a finite stopping power of the nu- 
clear matter. It simulates the finite stopping power by 
means of friction between three fluids. Naturally, this 
friction changes when deconfinement happens. In the 
case of the crossover scenario this change in the friction 
is very smooth. Therefore, it does not contribute to the 
irregularity. At the same time this change in the friction 



11 



enhances the irregularity in the first-order-transition sce- 
nario. As it was demonstrated in Ref. [ll], if the same 
friction is used in both phases, the reduced curvature cal- 
culated with the 2-phase EoS exhibits only a weak wiggle 
in Cy with considerably smaller amplitude as compared 
with zigzag in actual calculations with different frictions 
in different phases. These different frictions appear quite 
naturally in the 3FD model. The hadronic friction was 
estimated in Ref. [23| and works well at lower AGS en- 
ergies. Therefore, there are no reasons to modify it. The 
partonic friction, while not microscopically estimated, is 
fitted to reproduce data at high incident energies. This 
is a reason to believe that it is a proper choice. 

It is important to emphasize that the "peak-dip-peak- 
dip" irregularity is a signal from the hot and dense stage 
of the nuclear collision, rather than from the freeze-out 
stage as the most part of signal arc. 



VI. SUMMARY 

Proton and net-proton rapidity distributions in colli- 
sions of heavy nuclei Au-|-Au (at AGS and RHIC ener- 
gies) and Pb-|-Pb (at SPS energies) were analyzed in a 
wide range of incident energies 2.7 GeV < y'sjvw < 39 
GcV in terms of the center-of-mass energy per nucleon 
pair. The analysis was done within a model of the three- 
fluid dynamics employing three different equations of 
state: a purely hadronic EoS and two versions of EoS 
involving the deconfinement transition [l^ . These are an 
EoS with the first-order phase transition and that with 
a smooth crossover transition. The crossover transition 
constructed in Ref. [l^l is very smooth. The hadronic 
fraction survives up to very high energy densities. In this 
respect, this version of the crossover EoS certainly con- 
tradicts results of the lattice QCD calculations, where a 
fast crossover, at least at zero chemical potential, was 
found. Therefore, a true EoS is somewhere in between 
the crossover and 2-phasc EoS's of Ref. [l^ . 

Scenarios based on EoS's with deconfinement transi- 
tions have a theoretical advantage as compared to the 
purely hadronic one. In order to reproduce the baryon 
stopping at high incident energies, the friction between 
counter-streaming fluids have to be enhanced within the 
hadronic scenario as compared to its estimate based on 
experimental inclusive proton-proton cross sections j23j . 
Though such enhancement is admissible in view of un- 
certainties of the estimated friction, the value of the en- 
hancement looks too high. In scenarios with dcconflnc- 
mcnt there is no need to modify the hadronic friction. 
This can be considered as an indirect argument in favor 
of such scenarios. 

It was found that predictions within the flrst-order- 
transition scenario, i.e. with the 2-phase EoS, exhibit a 
"peak-dip-peak-dip" irregularity in the incident energy 
dependence of the form of the net-proton rapidity distri- 
butions. At low energies, up to Eiab = GeV, rapid- 
ity distributions for central collisions have a peak at the 



midrapidity, similarly to results with other EoS's. Be- 
ginning from 8A GeV, this peak turns into a dip at the 
midrapidity. Then again a peak is realized, starting from 
30A GeV. With further energy rise {Eiab > 40A GeV) 
the midrapidity peak again transforms into a dip, which 
is already a normal behavior which takes place at arbi- 
trary high energies. This behavior is in contrast with 
that for the hadronic scenario, where the form of distri- 
bution in central collisions gradually evolves from peak at 
the midrspidity (at Eiab < WA GeV) to a dip (at Eiab > 
10 A GeV). The case of the crossover EoS is intermediate. 
Only a weak wiggle of the type of "peak-dip-peak-dip" is 
observed in the energy range of lOA GcV < Eiab < 20^ 
GeV. 

The behavior the type of "peak-dip-peak-dip" in cen- 
tral collisions within the 2-phase-EoS scenario is very ro- 
bust with respect to variation of the model parameters in 
a wide range. It certainly disagrees with data at 8A GeV, 
10 A GcV and 40 A GcV energies. It also disagrees with 
data at 20A GcV and 30^ GeV, which have preliminary 
status, and hence it is too early to draw any conclusions 
from comparison with them. 

However, the experimental data also exhibit a trend of 
the "peak-dip-peak-dip" irregularity in the energy range 
8A GeV < Eiab < 40A GcV. Again this trend is based on 
preliminary data at energies of 20A GeV and 30^ GeV. 
Therefore, updated experimental results at 20A and 30A 
GeV are badly needed to pin down the preferable EoS and 
to check the trend of the "peak-dip-peak-dip" behavior in 
net-proton rapidity distributions. An irregularity in the 
baryon stopping is a signal of deconflnement occurring 
in the compression stage of a nuclear collision. It is a 
combined effect of the softest point of a EoS and a change 
in the nonequilibrium regime from hadronic to partonic 
one. It is important to emphasize that this irregularity 
is a signal from the hot and dense stage of the nuclear 
collision. 

An effective method to quantify the "peak-dip-peak- 
dip" irregularity is the analysis of the distribution shape 
in terms of the reduced curvature of the spectrum in the 
midrapidity Cy. In tems of Cy this irregularity is dis- 
tinctly seen as a zigzag irregularity in the energy depen- 
dence of Cy. 

Comparison with available data, including those 
at RHIC energies, indicate certain preference of the 
crossover EoS, though the crossover-EoS scenario in the 
presently used version does not perfectly reproduce the 
data either. Predictions of different scenarios for net- 
protons diverge to the largest extent in the energy region 
8A GcV < Eiab < 40A GeV. The prehminary data at 
20A and 30^ GeV disagree with any of the considered 
scenarios. At incident energies above the top SPS one 
the hadronic scenario certainly overestimates the proton 
midrapidity density, while deconflnement scenarios rea- 
sonably reproduce it. 

Anticipating results of subsequent papers, it should be 
mentioned that the 3FD simulations with same set of pa- 
rameters described here also reproduce other observables. 
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FIG. 9: (Color online) Midrapidity densities of various pro- 
duced particles as functions of the incident center-of-mass en- 
ergy of colliding nuclei predicted by 3FD calculations with 
different EoS's. Experimental data are from compilation of 
Ref. [5C Icomplemented by recent data from STAR collab- 
oration 53] and latest update of the compilation of NA49 
numerical results 1591. [6011. 



of course, with different degree of success depending on 
applied EoS. As an example, Fig. [HI demonstrates excita- 
tion functions of midrapidity values of various produced 
particles. As before, results for top calculated energy 
of y^SNN — 62.4 GeV should be taken with care, since 
accurate computation is still unavailable for this energy. 
The strangeness production at low incident energies is 
overestimated within the 3FD model. This is not sur- 
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FIG. 10: (Color online) Strangeness suppression factor as 
function of the center-of-mass energy of colliding nuclei. 



prising, since any EoS in the 3FD model is based on the 
grand canonical ensemble. This shortcoming can be eas- 
ily cured by introduction of a phenomenological factor 75 
[6l|, which accounts for an additional strangeness sup- 
pression due to constraints of canonical ensemble. The 
midrapidity densities of strange particles, displayed in 
Fig. ini are multiplied by 7s factor, which in its turn is 
presented in Fig. [TUl As seen, at Ei^h > WA GeV there 
is no need for additional strangeness suppression. 

As seen from Fig. ^ the purely hadronic EoS certainly 
fails at high energies. A preferable EoS is the crossover 
one, similarly to that for proton rapidity distributions. 
This calculations, as well as those of other observables, 
will be discussed in subsequent papers. 
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